Long wavelength local density of states oscillations near graphene step edges 
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Using scanning tunneling microscopy and spectroscopy, we have studied the local density of states 
(LDOS) of graphene over step edges in boron nitride. Long wavelength oscillations in the LDOS 
are observed with maxima parallel to the step edge. Their wavelength and amplitude are controlled 
by the energy of the quasiparticles allowing a direct probe of the graphene dispersion relation. We 
also observe a faster decay of the LDOS oscillations away from the step edge than in conventional 
metals. This is due to the chiral nature of the Dirac fermions in graphene. 



Defects, impurities and edges have a large effect of the 
density of states in metals and semiconductors. Their 
spatial extent is usually limited by screening which sup- 
presses the long-range behavior of the potential. How- 
ever, at the fermi wavevector, kp, the electron gas has 
a singularity leading to strong and long range oscilla- 
tions in the electron density known as Fricdel oscilla- 
tions. Fricdel oscillations have been observed near step 
edges on the surface of many metals such as CuQ and 
Au[2|- These materials are described by Fermi liquid the- 
ory and the oscillations have a characteristic wavelength 
given by twice the Fermi wavevector and a slow power 
law decay with distance. Graphene is a two-dimensional 
material where the quasiparticles obey the Dirac equa- 
tion and therefore behave as massless Dirac fermions |3|. 
Many unusual effects ranging from a different quantiza- 
tion of the quantum Hall effect 0, [Hj to the suppression of 
backscatteringQ have been observed in graphene. These 
effects are due to the chiral nature of the quasiparticles in 
graphene as well as the linear dispersion relation. These 
two properties also affect the Friedel oscillations and are 
expected to cause a faster decay with distance @. 

There have been several scanning tunneling 
troscopy studies of graphene on SiC0,lll and Si02[ 
which have shown interference patterns due to quasipar- 
ticle scattering. These studies have observed standing 
wave patterns due to scattering from defects in both 
monolayer and bilayer graphene and verified that the chi- 
ral nature of the quasiparticles leads to the suppression 
of backscattering for monolayer graphene Q. There has 
also been a topographic study of the effect of an armchair 
edge on the density of states in graphene on SiC which 
shows an oscillation along the C-C bonds that docs not 
depend on energy 
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However, none of these studies has 
observed a long wavelength oscillation in the density of 
states controlled by the energy of the quasiparticles along 
with its spatial decay. In this paper, we use scanning tun- 
neling spectroscopy to directly probe the local density of 
states (LDOS) near step edges in monolayer graphene. 
We observe a long wavelength oscillation of the LDOS 
due to scattering from the step edge. The wavelength of 
the oscillation is controlled by the energy of the quasi- 



particles. We also observe a faster decay {y^'^^'^) of these 
oscillations as compared to noble metals (y"^^^) due to 
the chiral nature of the quasiparticles in graphene. 

The measurements were done in an ultrahigh vac- 
uum, low-temperature STM (Omicron LT-STM) that 
was cooled to 4.5 K. In order to observe the long wave- 
length standing waves in the density of states, a very 
flat surface with few scattering centers other than the 
step edges was needed. Exfoliated graphene on Si02 
has about a 1 nm rou ghn ess due to the roughness of 
the underlying oxide 12, ll3[ . In order to have exfoliated 
graphene laying flat on a substrate we used hexagonal 
boron nitride (hBN) as a substrate. The sample was fab- 
ricated in a similar manner to previous graphene on hBN 
work 14 , 31 ■ Step edges in hBN occur naturally during 
the sample preparation procedure. The graphene lies flat 
on the hBN causing steps in the hBN to be transferred 
to the graphene. 

Figure [ija) shows a schematic of the experimental 
setup. The graphene flake is held at ground and the 
voltage on the tip is changed with respect to it. The 
tunnel current is measured from the tip to the sample. 
Figure [lljb) shows a 100 x 100 nm^ topography image 
with a step edge running almost horizontally. The step 
height is about 0.5 nm. On the higher portion (upper 
portion of Fig. [ijb)) of the step edge, a 10 nm Moire 
pattern can be seen. This is due to the lattice mismatch 
and rotation between the graphene and the underlying 
hBN as was reported earlier [ll,[3l. Figure [Ijc) shows a 
zoom-in on the lower portion (bottom of Fig. [Ub)) also 
showing a Moire pattern, but with a shorter ~2 nm pe- 
riod. The different Moire patterns are a clear signature 
of a relative rotation between the graphene and hBN lat- 
tices. An atomically resolved image (Fig. [TJd)) across 
the step edge shows that the graphene is continuous and 
therefore the step edge is formed in the hBN layer. The 
0.5 nm step is probably caused by an extra piece of hBN 
with a different lattice orientation. This different orien- 
tation leads to the two different Moire patterns on either 
side of the step edge. 

To study the influence of the step edge on the electronic 
properties of graphene, we have performed scanning tun- 
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FIG. 1. (color online) (a) Schematic illustration of the mea- 
surement setup showing the Au electrodes and graphene flake 
on top of a stepped hBN flake, (b) 100 nm x fOO nm to- 
pographic image of a step edge in the boron nitride beneath 
the graphene. (c) Atomic resolution image of the graphene 
lattice on the bottom of (b) showing a ~2 nm Moire pattern 
due to the underlying hBN. (d) Atomic resolution image of 
the graphene going over the step edge showing that it is con- 
tinuous. All of the topography images were acquired with a 
tip voltage of -0.3 V and a tunnel current of 100 pA. 



neling spectroscopy (STS) measurements over the area 
near the step edge. In these measurements, we recorded 
dl/dV; which is proportional to the LDOS, as a function 
of position and energy. The tip height was first stabi- 
hzed with a tip voltage of -0.3 V and 100 pA. Then the 
feedback circuit was turned off and a 5 mV ac voltage at 
563 Hz was added to the tip voltage. dl/dV as a func- 
tion of tip voltage was recored using lockin detection in 
order to create images of the LDOS. Figures ma)-(f) are 
maps of dl/dV for different tip voltages acquired in the 
region shown in figure [TJb). Two interesting features are 
noticed. First, LDOS oscillations in the lower portion 
arc present and they decay moving away from the step 
edge. There are also LDOS oscillations in the upper por- 
tion but the signal is significantly weaker. Second, the 
wavelength of the oscillations increases as the tip voltage 
decreases approaching the Fermi energy and Dirac point. 

LDOS oscillations have been previously observed on 
Cu[H and AuH surfaces. However, there is an obvious 
difference from our case. The wavelength of standing 
waves on noble metals is an order of magnitude shorter 
compared with that of graphene (1 nm compared to 10 
nm). Wc will see later that this difference comes from the 



different band structures. Long wavelength LDOS oscil- 
lations have been observed in other systems where the 
Fermi surface is not near the F point [l3, However, 
the chiral nature of the quasiparticles and suppression 
of backscattering in graphene leads to a different decay 
behavior of the oscillations as compared to these other 
materials. Recently, similar effects were also reported on 
topological insulators [l^ and explained by theoreti- 
cal calculations [2lj. Since both graphene and topological 
insulators are Dirac systems, they share a similar linear 
dispersion relation and chirality of electrons. 

We adapt the theory for topological insulators [2ll to 
our case of graphene. A step edge in graphene due to cur- 
vature can produce scalar and vector potentials which 
act as scatterers for electron waves 221. We define the 



direction along the step edge as the x axis (See figure 
[Sfa)). Consider an incident electron with wavevector 
K — Koirac + k whcrc Koirac IS thc Crystal momentum 
at the Dirac point and \k\ = fc is the radius of the con- 
stant energy circle in the Dirac cone. This electron will 
be partially reflected to a wavevector K' = K Dirac + k'. 
The electron will remain in the same Dirac cone as long 
as the step edge is not atomically sharp (the height of the 
graphene in Fig. [TJd) changes over a distance of about 2 
nm) and intervalley scattering is forbidden. A scattering 
barrier running in the x direction conserves the x com- 
ponent of k in the scattering process. So ii k ~ + ky 
then k' = k^ — ky in order to conserve energy and mo- 
mentum (see Fig. [3] (a) for illustration). The incident 
and scattered wave interfere with each other, and give 
rise to a spatially modulated electron density. So far we 
have only considered a specific pair of wavevectors, k and 
k' . In the actual scattering process, there are a range of 
wavevectors which are possible. All of the possible in- 
coming and scattering wavevectors k lie on a constant 
energy circle (CEC) determined by the tip voltage. In 
order to find the total change in the electron density, we 
must integrate the contributions from all of them, which 
gives. 



where ^{kx,ky) 



p(E,y)=f \^,{k,,ky)\^dk 

J CEC 



flection coefficient. Since the electrons in graphene are 
Dirac fermions they obey the 2D Dirac equation and the 
wave function in momentum space near the Dirac point 
is given by 



where 9k = arctan(-|ai). After reflecting off the barrier, 

Ky 

the phase factor changes io 6'^ ~ n — Ok- Then the spa- 
tially varying portion of the density of states is given by 

4r 

1^ cos{2kyy) sm{6k)dk 
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FIG. 2. (color online) (a)-(f) Images of the local density of states in graphene across a step edge running horizontally as a 
function of tip voltage. The tip voltages are (a) 138 mV (b) 118 mV (c) 98 mV (d) 78 mV (e) 58 mV and (f) 38 mV. The scale 
bar in all of the images is 10 nm. For all of the images, the tip height was stabilized at a tip voltage of -0.3 V and a current of 
100 pA. 



For large distances, the integrand varies rapidly and the 
only points that contribute are when ky is approximately 
constant. This occurs if ky is near its maximum value 
and hence the scattering is at nearly normal incidence 
(see Fig. [3] (a)). If the step edge is modeled as a po- 
tential barrier, then the reflection coefBcient is given by 
r oc sin(6'A:) for small Qu 2l|. Using the method of sta- 
tionary phase, the integral can be evaluated. Making the 
substitution, E = tmpk where Wi? is the Fermi velocity, 
we obtain an expression for the variation in the LDOS as 



6p{E,y) (X 



37r/4) 



(iJy)3/2 



(1) 



This expression describes the LDOS as a decaying oscil- 
lation with minima spaced by irhvp/E. Since we only 
probed energies near the Dirac point, which means the 
radius of the CEC and hence E is small, we get long 
wavelength LDOS oscillations compared with that from 
noble metals, where E is usually much larger. Another 
interesting feature of this expression is the decay rate. 
In noble metals, the LDOS oscillations decay with dis- 
tance from the step edge as I/y^Q. However, due to 
the pseudo-spin degree of freedom in graphene, backscat- 
tering is forbidden. This gives a faster, decay rate 

to the LDOS oscillation. 

Now we compare the experimental data with the theo- 
retical prediction. Figure |3](b) plots dl/dV as a function 
of tip voltage and distance away from a step edge. The 
data is acquired by averaging, in the x direction, all of the 
data in the two-dimensional images in figure [2l leaving a 
one dimensional plot of the LDOS versus distance. There 
has been an overall linear slope removed in the dl/dV 
curves as a function of energy. The data then shows os- 
cillations in the density of states with blues areas having 
higher density of states and red areas being low. The 



green and yellow lines show the minima in the density 
of states oscillations as predicated by equation ([1]). The 
only free parameter is the value of the Fermi velocity. 
We find slightly different values on the bottom (negative 
distances) and top (positive distances) of the step. For 
the yellow lines we found vp « 0.55 x 10^ m/s and for the 
green lines we found vp ~ 0.75 x 10^ m/s. Figure [3] (c) 
shows similar results from a different location which gives 
Vp ~ 0.45x 10^ m/s for the yellow lines and vp sa 0.5x 10^ 
m/s for the green lines. These values are lower than other 
measurements of the fermi velocity in graphene by about 
a factor of twoj^ [2^ [2^. The low values are confirmed 
by gate dependent spectroscopy measurements near the 
step edge which also give a reduced value of the Fermi 
velocity. The reduced value of the Fermi velocity may be 
due to the Moire pattern, as theoretical predictions have 
indicated that the fermi velocity can be strongly renor- 
malized for 2D periodic potentials which cause the Moire 
pattern As we observe two different Moire patterns 
on top of the step and below it, this may explain the dif- 
ferences in Fermi velocity observed. The reduced Fermi 
velocity may also be due to strain in the graphene caused 
by the step edge [2^. 

In addition to the oscillatory nature of the LDOS, we 
also expect it to decay with distance from the step edge as 
l/?/3/2. Figure Eld) plots dl/dV as a function of distance 
for a series of 5 tip voltages ranging from -30 mV to -90 
mV in steps of -15 mV. These curves also show fits (in 
blue) of equation ([T]). The curves show the decay 
expected for graphene in contrast to the l/y^^^ seen in 
conventional metals. From these fits we can extract both 
the amplitude and wavevector of the oscillations. 

Figure IDJa) plots the measured wavevector as a func- 
tion of tip voltage for the data in figure|3Ic). The dl/dV 
curve at each energy was fit using equation ([T]) to obtain 
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FIG. 3. (color online) (a) Schematic diagram of the scatter- 
ing process. Quasiparticles near the Dirac point are scattered 
from k to k'. In this process the x component of the wavevec- 
tor must be conserved because the step edge is along the x 
direction, (b) Density of states as a function of tip voltage 
and position. The data is averaged over the images in figure 
(2] Blue areas are high density of states while red represents 
low density of states. The colored lines represent the expected 
minima in the density of states from the theory calculations. 

(c) Density of states as a function of tip voltage and position 
for a difi'erent location showing similar behavior as in (b). 

(d) dl/dV as a function of position for several different tip 
voltages showing the decay with distance. The black curves 
are the experimental data and the blue curves are from the 
theoretical calculations. 
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FIG. 4. (a) Wavevector versus tip voltage extracted from 
fitting the dl/dV curves in figure [3jc). The solid line is a 
linear fit as expected for the graphene dispersion relation, 
(b) Amplitude of the LDOS oscillations as a function of tip 
voltage (energy). The solid lines are as expected from 

theory. 



consequence of the chiral nature of the quasiparticles in 
graphene. 
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the wavevector and amplitude at the given energy. We 
see that the wavevector is linear in tip voltage as expected 
from the dispersion relation of graphene giving a value of 
vf = 0.50 ± 0.05 X 10^ m/s. This linear relationship jus- 
tifies the fits in figure |3l^b) and (c) that assumed a linear 
relationship between E and k and used the fermi velocity 
as the fitting parameter. We can also plot the amplitude 
of the oscillations as a function of tip voltage which is 
shown in figure Ufb). From equation ([1]) we see that the 
amplitude of the oscillations is expected to decrease with 
increasing wavevector and hence energy. This is what is 
observed in figure |31[b) which shows a decrease of ampli- 
tude at larger tip voltages. The solid line is the expected 
decay of the amplitude as 1/ E^^'^ which agrees well with 
the measured data. 

In conclusion, we have demonstrated LDOS oscilla- 
tions near step edges in graphene on hBN. These os- 
cillations are long wavelength and are controlled by the 
energy of the quasiparticles. The decay of the oscilla- 
tions is much faster than in fermi liquids and is a direct 
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